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Abstract 

Bound states generated by the K coupled PT— symmetric square wells are studied 
in a series of models where the Hamiltonians are assumed 7?-— pseudo-Hermitian and 
7?.^— symmetric. Specific rotation-like generalized parities TZ are considered such 
that TZ^ = / at some integers A^. We show that and how our assumptions make the 
models exactly solvable and quasi-Hermitian. This means that they possess the real 
spectra as well as the standard probabilistic interpretation. 
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1 Introduction 



Bender's and Boettcher's PT— symmetric version of Quantum Mechanics [1] admits 
a transition to complex potentials V{x) (say, on a finite interval) characterized by 
the T'T— symmetry property VTV{x) — V{x)VT where V denotes parity while 
the complex conjugation T mimics time reversal. One of the simplest illustrative 
examples of the corresponding non-Hermitian PT— symmetric oscillator with real 
spectrum is generated by the purely imaginary square-well potential step [2] 

V{x) ^V^z){x) ^-iZsign{x), x e (-1, 1). (1) 

The solvability of this model facilitates an introduction of the norm of its wave func- 
tions in a suitable metric, i.e., via an introduction of a Hamiltonian-dcpendent scalar 
product [3]. It also renders possible the correct transition to classical limit [4]. In- 
teresting physical applications of eq. (1) were found in supersymmetric context [5] 
as well as beyond quantum theory where schematic eq. (1) and its modifications 
may play role in an explanation of the mode-swapping phenomena in classical mag- 
netohydrodynamics [6]. In mathematical context, model (1) helped to clarify the 
mechanisms of the spontaneous breakdown of PT— symmetry at a critical strength 
Z — Zcrit of non-Hermiticity [7, 8]. 

A specific merit of model (1) lies in the feasibility of a transition to its more 
sophisticated piece- wise constant solvable alternatives [9]. A particularly promising 
new direction of development has recently been found in the tentative use of the ele- 
mentary functions (1) as forces which mediate an interaction between two [10] and/or 
three [11] coupled square- well oscillators. Here we intend to move one step further 
and to re-analyze the similar coupled channel problems in a more systematic 
manner. 

Our key idea lies in the observation that in the one-dimensional Schrodinger 
equation which describes K coupled channels, 

"^^^"^^"^^ + t{ ^^(-^^ '^^'^^^^ ^ V-^H^), m = 1, 2, . . . , X, (2) 

the complications connected with its solution grow very quickly with K. In general, 
the properties of the model are controlled by as many as independent real cou- 
plings ^(mj) in eq. (1). An introduction of some additional symmetries would be 
desirable. Recently we successfully reduced the number of free parameters to three 
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in the K — 2 model of ref. [10] and, under the "stronger" symmetry assumptions, to 
two in the K — 3 model of ref. [11]. 

Inspired by the latter two examples we shall now contemplate K > 3 channels 
and try to impose certain symmetry constraints in the manner which could keep 
our Schrodinger equation (2) with more coupled channels exactly and compactly 
solvable. 

2 PT— symmetry revisited 

2.1 Parity re-interpreted as a pseudo- metric 

In the majority of its updated formulations [12, 13, 14], PT— symmetric Quantum 
Mechanics (PTSQM) replaces the involutive parity V — by an invertible and 
indefinite pseudo-metric P. In the language coined by Ali Mostafazadeh [13] one 
replaces the PT— symmetry property of the Hamiltonian H ^ by the requirement 

i/t = Pi/p-\ P = pt^7. (3) 

It may be understood as a certain necessary condition that the spectra of the ob- 
servable H and of its "redundant" conjugate H'^ coincide. 

A "hidden" purpose of the postulate (3) lies in the fact that as long as H ^ H\ 
the standard knowledge of the (presumably, real and discrete) energies En and of the 
related eigcnstatcs \n) of H must be complemented by the independent construction 
of the eigenstatcs of the conjugate operator H\ i.e., in our adapted Dirac's notation, 
of the "ketkets" \n)) . Fortunately, in such a situation eq. (3) enables us to employ, 
in the non-degenerate case, the implication 

H'^ \n)) = En \n)) =^ |n)) = const{n) P \n) . (4) 

The latter condition is of paramount importance for the technical feasibility of the 
practical applications of the formalism. Indeed, the construction of \n)) becomes 
straightforward whenever the action of the pseudo-metric P is not too complicated. 

A more detailed support of the latter argument may be found, e.g., in the Appen- 
dices of ref. [10] and in ref. [15] where we studied an apphcation of PTSQM to the 
Peano-Baker-like two-channel version of Klein-Gordon equation. On the background 
of this technical summary we only have to emphasize that the knowledge of the two 
sets of the vectors \n) and \n)) (forming a biorthogonal basis in our Hilbert space) 
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opens, for all the square-well-type models, the way towards the construction of the 
necessary "physical" positive definite metric © = A nice explicit illustration of 
the recipe (which ascribes an appropriate probabilistic interpretation to the system, 
cf. Appendix A below) has been discussed by Mostafazadeh and Batal [4], with an 
elegant Krein-space mathematical re-interpretation added recently by Langer and 
Tretter [8]. 



2.2 Unitary alternatives to the parity 



We feel inspired by the observation that in a perceivable contrast to the physical 
metric © itself, the pseudo-metric plays just an auxiliary role, via eq. (4). In such a 
context, the requirement of the Hermiticity of P is redundant and the main emphasis 
must be put on its simplicity. This is the key idea of our present paper. In place of the 
standard Hermitian pseudo-metrics P we shall try to work with some non-Hermitian 
pseudo-parity operators R 7^ replacing eq. (3) by its alternative 



(5) 



As long as the Hermitian conjugation is an involution, we may insert cq. (5) in its 

conjugate version H = [R~^]^ H'^ and arrive at the symmetry requirement 



HS^SH, 5 = [R-^]^R. 



(6) 



We believe that its implementation might serve our present purposes. 

For the sake of definiteness we shall pay attention to the families of operators R, 

^ ... v\ f ... V ^ 



R, 



iK,l) 



V 
V 

















R 



{K,2) = 




V 












V 

V 

... 





yO...O V J yu... u uy 

with the property R~^ = R^ giving the symmetry S — R^. In the representation 
where we separate the parity, R{k,l) — 'Py(^k,l)-i we may employ the recurrences 



r(x,L+i) = r(x,i) Y^K,L) , 
At any K the pseudoparities obey the rule 



"(X,L) 



K 



L = l,2,... . 

= / and may be interpreted as 

t 



finite rotations. At the even K — 2M we get the non-Hermitian yi^k,k-l) 



■iK,L) 
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ai L — 1, 2, . . . , M — 1. An anomaly occurs ai L — M where we note that r(^K,K-L) — 
^{K,L) ^ — r(2M,M) remains Hermitian. There is no similar Hermitian exception at 
the odd integers X = 2M - 1 with M > 1. 

We are now prepared to apply the rule (5) to our coupled-channel models (2). 



3 The method 

3.1 The symmetry-compatible sets of coupling constants 

In units H — 2m — 1 and in a partitioned matrix notation equation (2) may be 
reformulated as a diagonalization of the Hamiltonians H — H(^kinetic) + H (interaction) {x) 
where 

/ _^ n ^ 



■-. 



(kinetic) 







V 











- 



The separate parity-antisymmetric imaginary square- well potentials (1) will form the 
array 



H (inter action) (•^) 



V ^Z^K,l)i^) ^2(K,2)(^) 

characterized by K"^ different real couphng constants, 



(K,K) 



[x) 



(7) 



^(1,1) ^(1,2) 



Z{2,\) ^(2,2) 



Z(l,K) 
Z{2,K) 



\ Z(K,1) Z(^K,2) ■ ■ ■ Z(^K,K) ) 



(8) 



As long as the diagonal H^j^inetic) commutes with all our pseudoparities, we shall 
only have to study the consequences of the symmetry (6) upon the variability of the 
set (8). 

In a more detailed analysis of the latter point we must recollect, firstly, the action 
of parity V on our Hamiltonian, 



VH, 



(interaction) 



(interaction) 
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This rule enables us to re-write eq. (5) as a condition imposed upon the real matrix 
of indices (8), 

A = r{K,K-L) ■ • ri^K,L) (9) 

where denotes transposition. This equation is a core of our forthcoming construc- 
tions. At any fixed K and L, it must be satisfied as a guarantee that our square-well 
Hamiltonian H obeys the R— pseudo-Hermiticity rule (5). Step-by-step we shall list 
the solutions of eq. (9) distinguishing between the odd K (in a series starting at 
section 4) and even K (starting from section 6). 

3.2 The determination of the bound-state energies 

Our next step may be guided by the elementary single-channel example of ref. [2] 
with the "effective" Schrodinger equation 

--^^{x) -^Zeffsign{x)ip{x) = Eip{x), (^(-1) = (^(1) = . (10) 

At K — 1 and Zeff — Z it has been shown solvable and physical (i.e., possessing 
the real spectrum) at |^e//| < -^cr-it ~ 4.48 in ref. [2], with a better estimate of 
Z^rit ~ 4.47530856 derived in ref. [11]. 

For the generic K > 1 and for the negative or positive coordinate x, our coupled 
set (2) are differential equations with constant coefficients. In a way resembling 
eq. (10) these equations remain solvable by the trigonometric ansatz 

1 C'f^ sinfi:^(a; + l), xG (-1,0), 

\ siuKRi-x + l), a;e(0,l), m = l,2,...,X 

compatible with the "external" boundary conditions at x — ±1. We must also 
impose the 2K—p[et of the standard "internal" matching conditions in the origin, 

C^^^ sin kl = Clr^ siuKR, m = 1, 2, . . . , X, ^^^^ 

kl C^l^ cos kl = -kr C^r^ coskr, n^l,2,...,K. 

Their first half may be read as determining, say, the K "dependent" constants Cj^^ 
as functions of the K "independent" constants C}"*'' and of the not yet specified 
parameters k^^r. The ratio of the equations with m = n eliminates all the constants 
and leads to the single complex condition 

KLCOtSiO-KL — —KRCOtanKR. (13) 
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Finally, the insertion of our ansatz (11) in the Schrodinger equation (2) gives 



1^1 + i ^(1,1) i Z(i,2) 
i ^(2,1) 



4 + iz, 



(2,2) 



iZ, 

i^i 



\ 



V 



/ ..2 



c 



E 



at a; G (—1,0) and the similar i^— dimensional diagonalization 



{K) 



(14) 



-iZ 



(2,1) 



— iZ| 



(2,2) 



-iZ 



(2,i^) 



v 



-iz. 



-iZ, 



(K,2) 



K 



R 



C 



(2) 
R 



\^R J 



c 



E 



c 



R 

(2) 
R 



\^R J 



(15) 



'{K,K) j 

at X G (0, 1). For the real energies the latter two sets are just complex conjugates of 
each other so that we may omit one of them and set — s ^ it — k\ with s > 
and t G (— oo, oo). 

In the next step generalizing the experience gained in refs. [2, 10, 11] we use 
another ansatz 

E^s^ - 1^. (16) 

Its use leaves the matrices in eqs. (14) and (15) purely imaginary so that it is easy 
to write down their common secular equation which is real. 



det 



— '2-St -^(1,2) 

-Z'(2,i) -^(2,2) — 2si 



Z(1,K) 
Z{2,K) 



\ 



0. 



(17) 



y Z{K,l) ^{K,2) ■ ■ ■ Z(^K,K) — 'ist j 

It is to be complemented by the complex eq. (13) which, in the same notation, 
degenerates to the semi-trigonometric and K— independent algebraic formula 



2s sin 2s + 2t sinh 2t = 0. 



(18) 



We may summarize that the pair of the real algebraic equations (17) and (18) may 
be expected to specify the real parameters s = s„, t — tn and the energies E — Em 
n = 0, 1, ... at each particular number of channels K. 

As long as equation (18) itself is the same for any we may treat it simply as a 
definition of certain "universal" curve t = texacti^)- It is worth noting that its shape 



6 



carries a lot of resemblance to the half-ovals 

tauxiiiary{s) = max(0, s sin 2s) . 

(see ref . [2] for a more detailed description of the shape of the exact curve) . 

The former equation (17) defines, in principle, the dependent K—plet of the 
real eigenvalues 2st = Zgff. Thus, once we determine all of them exactly, Zgff — Z^ff, 
k — 1,2, . . . , K, we may interpret all the related physical roots s — Sn and t — tn 
(which may be degenerate of course) as the coordinates of the intersections of the 
above-mentioned universal half-oval curve t — texact{s) with any one of the K much 
more elementary hyperbolic curves 4wboiic(«) = ^e///(2s)- 



4 The first nontrivial model with odd = 3 



4.1 Hermit ian choices of pseudo-metrics P 

When we pick up X = 3 we have a nice opportunity to distinguish between the 
parity V, Hermitian pseudometric P and its non-Hermitian generalization R. In the 
simplest Hermitian arrangement we may choose the diagonal partitioned operator 

^ P \ 
P = P(3,o)= P = Pt = p-i 
^0 Q V J 

and notice that its use in eq. (5) does not impose any constraint upon the nine 
one-parametric square-well interactions 

iZdk^. X e (-1,0), 

(0,1), 

All the nine coupling constants Z = Z^i^-^ remain independent. This leaves the 
corresponding solutions very complicated. There are no symmetries in the problem, it 
must be solved more or less purely numerically. This is, from our present constructive 
point of view, a not too interesting situation. 

The situation is merely marginally improved by the transition to several other, 
less trivial Hermitian pseudo-metrics like 

/ ^ 
P = P(3,i) = P = Pt = p-i 



'{3,k)-, 



j. A; = 1,2, 3. 



(19) 
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or 



(3,2) 



/ 7^ ^ 

V 
\0 V ^ 



It is easy to show that in the Hermitian cases the partitioned form of the second 



powers 



1 2 



(3,1) 



(3,2) 



IV^ becomes diagonal. 



4.2 Non-Hermitian pseudo-parities R 

The emerging availabihty of the first two simplest non-Hermitian pseudo-parities 
should be emphasized at K — 3, 



R = R 



(3,1) 



f v\ 

POO 
y V J 



R 



(3,2)' 



R — R 



(3,2) 



f V \ 
OOP 
POO 



^(3,1)' 



^ P \ 
R"^ = OOP 

^ P y 

/ P ^ 
R~^ = POO = R, 
y P y 



(3,1)- 



In contrast to the Hermitian cases, we only obtain the higher-power parallel rules 



R 



(3,1) 



jp3 and 



R 



(3,2) 



IV since 



R 



(3,1) 



PR, 



(3,2) 



and 



R 



(3,2) 



P R(3,i) in both our genuine non-Hermitian samples. 



4.3 The allowed ^(j,A;) for the three coupled channels 

At the unique index L = 1 our R— pseudo-Hermiticity condition (9) acquires the 
linear algebraic form of a set of equations for the couplings Zq^^) = a{j, k), 



/ 1 \ 
1 

U 1 0/ 



^ a(l,l) a(l,2) a(l,3) ^ 

a(2,l) a(2,2) a(2, 3) 
^ a(3, 1) a(3,2) a(3, 3) y 



^ a(l,l) a(2,l) a(3,l) ^ 
= a(l,2) a(2,2) a(3, 2) 

^ a(l,3) a(2,3) a(3, 3) j 
Due to the non-Hermiticity of the pertaining matrix r, these relations represent a 
much more powerful constraint which leaves just the two coupling constants free and 



^001^ 
1 

U 1 0/ 
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independent. The equations preserve their form under the transposition so that the 
structure of our K — 2 square-well model remains independent of L, 

( Z X x\ 

-^(interaction) = X Z X , L = 1,2. (20) 

yX X Z J 

Now we may return to the Hermitian choices of P = P(3,i) or P = P(3,2) of subsection 
4.1 which, obviously, introduced much less symmetry. Indeed, the solution of the 
corresponding nine linear equations generates just the three nontrivial constaints so 
that as many as six coupling constants remain independently variable. The same 
tendency survives at the higher K. We may conjecture that the breakdown P — > R 
of the Hermiticity is connected, definitely, with an enhancement of the symmetry 
and with the simplicity of the models. 



4.4 Energy levels 

In the final step of the construction of the bound states at X = 3 we may follow 
either the general recipe of section 3.2 or the recent detailed presentation of the 
K — 3 solutions in ref. [11]. In essence, we have to connect the parameters s and t 
with the "effective charge'", 

2st^Z,ff, Z,ff^Z + F 

where, in the notation of eq. (20), the three eligible values of the shift F — Fj are to 
be sought as eigenvalues [Fi, F2, F3] = [2X, —X, —X] of the modified matrix (20), 

^ X X^ 
X X . 

^X X J 

The three respective eigenvectors may be found in ref. [11] - here our MAPLE 
software produced their following simpler alternative sample 

{1,1,1}, {-1,0,1}, {-1,1,0} 

which is still to be re-orthogonalized. 

In a way compatible with the results of ref. [11] we may summarize that our K — 3 
coupled bound states are determined by formulae (11) and (16). The parameters s 
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and t are fixed as intersections of the half-ovals (18) with one of the two available 
hyperbolic curves, 



t = &\s) 



2s 



'eff 



[a), a = l,2, Z,ff{l) = Z + 2X, Z,ff{2) = Z - X . {21) 



By construction, the second family of intersections represents the twice-degenerate 
levels. 



5 The next model with odd K = b 



There is no anomaly in the non-Hermiticity of r = r(5_i) with L = 1,2,3 and 4, 



r(5,l) = 



/ 1 \ 

1 

10 

10 

V 1 ; 



r(5,4) = 



/ 1 \ 

10 

1 

1 

V 1 ; 



^{interaction) 



(22) 



All the four different pseudo-Hermiticity conditions (9) with L = 1, 2, 3, 4 lead to the 
same three-parametric coupling-constant matrix 

f Z X D D X \ 

X Z X D D 

D X Z X D 

D D X Z X 

\X D D X Z ) 

Besides its exceptional eigenvector {1, 1, 1, 1, 1} pertaining to the obvious eigenvalue 
Fq — 2 D + 2 X , the reduced Z = form of this matrix is most easily shown to 
possesses the pair of the twice degenerate eigenvalues, 

^ 2 

with the two respective eigenvectors 



-D-X±V5{-D + X) 



and 



"1 D-X±^/5{-D + X) 1 D-X±V5{-D + X) ^ _ 1 
2 -X + D '~2 -X + D ' ' '~ J 



1 D - X ±y/5{-D + X) 1 D-X±^/5{-D + X) 

2 -X + D ,-1,0,1,-- -X + D 
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The next steps towards the next odd K — 7,9, . . . will be discussed in our concluding 
remarks. Now, for a more specific illustration of some technical subtleties let us 
return to the systems with the small even numbers of coupled channels. 

6 The simplest model with even K = 2 

It would be easy to relax the involution assumption = / as formally redundant 
and pedagogically partially misleading. Even without such a constraint we get just 
the most elementary Hermitian option at = 2, 

plus its fully off-diagonal alternative 

\V J J 

Only for the latter sample choice of the matrix r = r(2,i) the compactified version 
(9) of the pseudo-Hermiticity condition (5) acquires a nontrivial two-by-two matrix 
form, 

1 W a(l,l) a(l,2) \f a(l,l) a(2, 1)) U 1 ^ 

1 ) [ a(2,l) a(2,2) ) [ a(l,2) a(2,2) )[l )' 

With involutive V this model as well as the resulting set of the four equations has 
already been studied in ref. [10] and may be easily shown to degenerate to the single 
constraint a(l, 1) = a(2,2). Thus, the K = 2 version of our present square-well 
model possesses three free real parameters X, Y and Z, 

■^(interaction) 

As long as the pseudo-parity remains Hermitian, the model does not fit in the scope 
of our present paper. Still it exemplifies the general pattern of the construction since 
both the eigenvalues F — F± oi the modified Z — version of matrix (23) are easily 
found, F± = ±\/XY and also the determination of the two respective eigenvectors 
by our MAPLE program, 

|i,^x7r|,{i,-^/x7r| 

is easily verified by hand and tests the recipe. 




(23) 
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7 Four coupled channels 



Just a smaller representative sample is to be added at = 4, viz, the non-Hermitian 



R — R'(4,i) — 



/ 
V 



V 



\ 





p y 



and the "exceptional" Hermitian 

/ P \ 
COOP 
P 



R — R 



(4,2) 



R^ R-^ = 



\ p y 

It should be noticed that while we obtain a diagonal 





/ 


p-1 








\ 


-1 _ 


























p-1 






\v-^ 











/ 




' 





p-1 





\ 


1 _ 











p-1 






p-1 















V 


p-1 








/ 



R 



(4,2) 



= 7 7^^ in the second 



power of our Hermitian operator, and analogous non-Hermitian formula requires the 

l4 



use of the fourth power. 



R 



(4,1) 



The overall structure of the matrix of couplings A compatible with our require- 
ment of the "maximal" non-Hermitian symmetries (9) ceases to be unique at A' = 4. 
A 'priori^ this follows from the observation that besides the matrix-transposition 
mapping between the pseudo-Hermiticity constraints at L = 1 and L = 3, one also 
encounters the anomalous Hermitian problem at L = 2. In the latter case we have 
to solve the 16 linear algebraic equations 

/O 1 0\/a(l,l) a(l,2) a(l,3) a(l,4)^ 

a(2,l) a(2,2) a(2, 3) a(2,4) 

a(3, 1) a(3,2) a(3, 3) a(3,4) 

V a(4,l) a(4,2) a(4, 3) a(4,4) j 



1 
10 

1 y 



o\ 
1 



/ a(l,l) a(2,l) a(3,l) a(4, 1) \ 

a(l,2) a(2,2) a(3, 2) a(4, 2) 

a(l,3) a(2,3) a(3, 3) a(4, 3) 

V a(l,4) a(2,4) a(3,4) a(4,4) j 

which only impose the six constraints upon the 16 couplings. Prom the practical 
point of view, too many of them remain freely variable. Similar results are obtained 
also for the other Hermitian operators R. 



/ 1 



10 

V 1 y 
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In contrast, the parallel and transposition-related L — 1 and L — 3 non-Hermitian 
versions of eqs. (5) give the set of 16 equations 



and 



/ 1 \ 

10 

10 

V 1 / 

/ 1 \ 

10 

1 

V 1 y 



A = 



/ 1 \ 

10 

10 

\^ 1 y 

/ 1 \ 

10 

1 

V 1 y 



L = 3, 



respectively. Both of them give the same four-parametric set of coupling constants 



^{interaction) 



( z u 

L Z 

D U 

\ L D 



D U \ 

L D 

Z U 

L Z J 



L = l,3, 



(24) 



compatible with both our non-Hermitian pseudo-parities and with the four quadru- 
plets of the independent elements sitting on the main diagonal (Z), side diagonals 
(D) and in an upper square (U) and lower square (L). After an obvious permutation 
of the basis this matrix may be understood as a partitioned structure 



i (permuted) 
(interaction) 



z 


D 


U 


u 


D 


Z 


u 


u 


L 


L 


z 


D 


L 


L 


D 


Z ^ 



L = l,3, 



(25) 



which represents an partitioned generalization of the two-by-two model (23) above. 
Directly, this solution may be derived from the pseudo-parity 

/ 1 \ 
1 



„(permuted) 'p—l-£^(permuted) 



10 

V 1 y 



which is just a permuted version of the above non-Hermitian square root r(4 1) of the 
unit matrix. 
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Four eigenvalues Fj — [—D, —D, D + 2 VUL, D — 2 VUL] of the modified matrix 
(24) with Z — form the partially degenerate quadruplet. The respective eigevectors 
read 



{1, 0, -1, 0} , {0, 1, 0, -1} , |l, 1, ^ j , {l, 1, 



u 



u 



u 



The quadruplet of eigenvalues derived from the alternative four-by-four matrix (25) 
with Z — remains unchanged of course. Even the respective eigenvectors them- 
selves become merely predictably influenced by the underlying permutation, 



{0,0, 1,-1}, {-1,1, 0,0}, 1,1, 



U 



u 



1,1, 



u 



/UL 

IT 



8 Six coupled channels 

While the choice of the Hermitian pseudo-parity r(6,3) is not sufficiently restrictive 
and leaves 21 free parameters in the related six-by-six coupling-matrix A, much more 
symmetry (with just 7 free parameters) is induced by all the non-Hermitian r(6,L) 
with L 3. 

It is worth noting that different patterns are obtained at L = 1 or L = 5 (when 
the resulting real matrix A remains asymmetric) and at L = 2 or L = 4 (when 
the resulting real matrix A becomes symmetric). In the former suitable 
permutation of the matrix indices leads to a maximally compact picture. 



(permuted) 
■(6,1) 



1 

1 

1 

1 

1 

1 



A 



(permuted) 
(interaction) 



\ 



z 


Y 


G 


B 


F 


B 


X 


Z 


C 


F 


C 


G 


F 


B 


Z 


Y 


G 


B 


C 


G 


X 


Z 


C 


F 


G 


B 


F 


B 


Z 


Y 


C 


F 


C 


G 


X 


Z 



\ 



(26) 

The partitioning of the latter asymmetric matrix indicates how our system may be 
visualized as a coupled set of its three two-dimensional asymmetric subsystems of 
the form (23). 

For the second option with L = 2 it is remarkable to notice that while the former 
pseudo-parity remains asymmetric and, hence, non-Hermitian, the latter matrix of 
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coupling constants A appears to be, for some unknown reason, symmetric. 



/ 



{permuted) 
'(6,2) 



1 

1 
1 
1 
1 
1 



(permuted) 
(interaction) 



z 


X 


X 


G 


D 


G 


X 


z 


X 


G 


G 


D 


X 


X 


z 


D 


G 


G 


c 


G 


D 


A 


B 


B 


D 


G 


G 


B 


A 


B 


G 


D 


G 


B 


B 


A 



(27) 

Moreover, its inspection reveals that it is again tractable as a coupled system of its 
two three-dimensional subsystems of the symmetric-matrix form (20). 

One has to add that the diagonalization of the Z — version of the matrix (26) 
is still feasible and gives the two nondegenerate eigenvalues 



^±0 



G + R±V2GY + 2BX + XY + ABG 



and the two doubly degenerate eigenvalues 



1 



'±1 



— G--R±- V-3G^ + 6GR-3R'^ -ABX + 4:XY + ABG -AGY. 
2 2 2 



In contrast, the study of the model (27) is hindered by the occurrence of the two 
different couplings on the main diagonal. One must employ a shift of the eigenvalues 
F — X + {Z + A)/2 which leads to the diagonalization of the matrix of the form 



V 



—00 


X 


X 


G 


D 


G 


X 


—CO 


X 


G 


G 


D 


X 


X 


—LU 


D 


G 


G 


G 


G 


D 


CO 


Y 


Y 


D 


G 


G 


Y 


CO 


Y 


G 


D 


G 


Y 


Y 





\ 



J 



with non- vanishing main diagonal, cu — {Z — A)/2. Still, in a way resembling the 
previous L — 1 model we get the two nondegenerate eigenvalues 



A 



±0 



Y + X±JAo 



with the abbreviation 



Ao = Y^ -2XY + X^ + co^+ 
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+2 ya; - 2 Xu; + + 2 DC + + 2 + + 2 GC, 
accompanied by the two doubly degenerate eigenvalues 

A±i 

where the discriminant Ai is equal to the sum 



2 2 2 ^' 



Y^-2XY + X^ + 4:D^ + 4:CU^-4: Yu+ 



+AXu;-4:DC + 4:C'^ -AGD-AGC + AG^. 



9 Concluding remarks 

9.1 Seven and more coupled channels at odd K = 2M — 1 

The pattern initiated by X = 3 and X = 5 is perpetuated at the next odd dimension 
K — 7 where the linear system of 49 equations yields the four free parameters in the 
L— independent solution 

/ y V V n n V v \ 



{interaction) 



V 



z 


X 


Y 


D 


D 


Y 


X 


X 


Z 


X 


Y 


D 


D 


Y 


Y 


X 


Z 


X 


Y 


D 


D 


D 


Y 


X 


Z 


X 


Y 


D 


D 


D 


Y 


X 


Z 


X 


Y 


Y 


D 


D 


Y 


X 


Z 


X 


X 


Y 


D 


D 


Y 


X 


Z 



(28) 



It is easy to sec the general pattern and to guess the structure of A (^interaction) at 
any higher odd dimension K — 2M — 1. At any M this hypothesis may be verified, 
say, in MAPLE, by using the same algorithm as employed in the above calculations 
performed at the first three nontrivial indices M — 2, M — 3 and M — A. 

In a way extrapolating the K — b results we may move to X = 7 in eq. (28) and 
guess the exceptional eigenvector {1,1,1,1,1,1,1} and its eigenvalue Fq — 2D + 
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2X -\-2Y . The three further eigenvalues of the matrix 



n 
u 


Y 


1 






1 


Y 


X 





X 


Y 


D 


D 


Y 


Y 


X 





X 


Y 


D 


D 


D 


Y 


X 





X 


Y 


D 


D 


D 


Y 


X 





X 


Y 


Y 


D 


D 


Y 


X 





X 


X 


Y 


D 


D 


Y 


X 






are doubly degenerate and may be written in the form of Cardano formulae. They 
represent the real roots in terms of the complex quantities and, for this reason, 
we omit them here and leave their generation to the interested readers as an easy 
exercise. 



9.2 Eight and more coupled channels at even K = 2M 

One gets lost when solving eq. (9) in the Hermitian case at L = 4 with 37 free 
parameters in A8 4, 29 of which occur there in pairs. 

Similarly, not enough symmetry is induced by the L — 2 pseudo-parity since the 
resulting asymmetric A depends on as many as 16 free parameters, each occurring 
strictly four times. 

Thus, the only satisfactory reduction is obtained from the non-Hcrmitian pseudo- 
parity r(8^L) with L = 1 (and, identically, L = 3 etc). The explicit solution of the 
corresponding eq. (9) leads to the permission of the eight free coupling strengths, 
each occurring eight times in A. One parameter sits simply on the main diagonal and 
one on the two "submain" diagonals of the two off-diagonal quadrants. The remaining 
six octets form the asymmetric pattern of three doublets reflected pairwise by the 
main diagonal. The first pair jumps over the main diagonal and the second one over 
the two submain diagonals while the last pair simply fills the remaining vacancies. 
We may again separate the even and odd rows and columns and obtain a re-arranged 
matrix A which is partitioned in the four submatrices where the two diagonal ones 
are composed of the four parameters only and represent just a transposition of each 
other. Each of the two off-diagonal blocks depends just on the two parameters in a 
way resembling slightly the partitioned structure of eq. (25). 
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Appendix A: A note on the interpretation of the 
Hamiltonians 



In a less usual approach to Quantum Mechanics as outlined in review [3] the Hamilto- 
nian is allowed to be quasi-Hermitian, i.e., Hermitian with respect to some "anoma- 
lous", nontrivial metric 7^ / in the Hilbert space of states. Of course, one must, 
first of all, re-define the new, "anomalous" scalar product 

\a) © \b) = {a\e\b) , e = > 0. (29) 

In effect, this is equivalent to a replacement of the standard Hilbert space H by its 
new and ©—dependent "physical" version V equipped with the more flexible and 
adaptable product (29). 

The key point is that once we stay within the innovated space V, all the basic 
principles of Quantum Mechanics remain unchanged. At the same time, the new 
flexibility carried by our freedom of the choice of © is compensated by the loss of 
the meaning of the standard Hermitian conjugation if — > if^ Indeed, the standard 
textbook Hermiticity of the observables A — must be replaced by the requirement 

\Aa) © \b) = \a) © \ Ab) 

which may be re-read as an obligatory quasi-Hermiticity property 

A^'^QAQ-^ (30) 

of all the observables A in the new formalism. Of course, no new physics is being 
discovered in this manner because all the operators A with the property (30) are 
simply Hermitian with respect to the fixed new metric © 7^ 7 [3] denoted usually as 
77+ by Mostafazadeh [13] and factorized as © = CP with "charge" C by Bender et 
al [14]. 
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